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The effect of nuclear deformation on the energy distributions of the Gamow- Teller strength is 

studied in neutron-deficient Hg, Pb, and Po even isotopes. The theoretical framework is based on 

a self-consistent deformed Skyrme Hartree-Fock mean field with pairing correlations between like 

nucleons in BCS approximation and residual spin-isospin interactions treated in the proton-neutron 

quasiparticle random phase approximation. After a systematic study of the Gamow- Teller strength 

\^ • distributions in the low excitation energy region, relevant for /3^-decay, we have identified the best 

f^ I candidates to look for deformation signatures in their /3"'"-decay patterns. /S"*" half-lives and total 

^^ ■ Gamow- Teller strengths B(GT^) are analyzed as well. 
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I. INTRODUCTION 



The existence of spherical and deformed nuclear shapes coexisting in a given nucleus at low excitation energies is 
I> I nowadays a well established feature characterizing many isotopes in the neutron-deficient Hg, Pb, and Po region (see 
\l Ref. [1] and references therein) . These coexisting structures have been found experimentally by studying the 7-rays in 

^^ coincidence with the emitted a-particles in the a-decay of the fusion products created in selected fusion-evaporation 
~l : reactions. As a matter of fact, the existence of at least one low-lying excited O''^ state in all even Pb isotopes between 
^-v ' A = 184 and A = 194 has been experimentally observed at excitation energies below 1 MeV [1,2]. Similarly, Hg 
\^ ' isotopes from A = 180 up to A = 190 have been found to be oblate in their ground states with prolate excited states 
(^ ■ in the 1 MeV range. 

^H ', Within a mean field description of the nuclear structure, the presence of several 0"'" states at low energies is under- 

-4— > . stood as due to the coexistence of different collective nuclear shapes. The energies of the different shape configurations 
,— H ' are calculated using a nuclear potential with the energies of the single particle orbitals depending on the deforma- 
y , tion. Calculations of the potential energy surface have become more and more refined with time. Phenomenological 
mean fields and Strutinsky method [3] , are already able to predict the existence of several competing minima in the 
deformation energy surface of neutron-deficient Hg and Pb isotopes. Self-consistent mean field calculations [4,5] and 
calculations including correlations beyond mean field [6-8] confirm these results. In particular, for Pb isotopes, all 
the approaches analyzed in Ref. [7], from mean field up to very sophisticated angular momentum projected generator 
coordinate methods, provide the same underlying basic picture of strong coexisting spherical and deformed shapes. 
C^ This justifies the use of mean field approaches as a first estimate to a qualitative description of the energy minima 

and to associate the 0+ states with coexisting energy minima in this energy surface. 

The relative energies of the different shapes predicted by mean field calculations in the neutron-deficient Hg-Pb 
region have been found to be very sensitive to fine details in the calculations, specially to pairing effects. This is 
true in both non-relativistic [9] and relativistic calculations [5,10,11]. In particular, standard deformed relativistic 
mean field (RMF) calculations [10,12] do not reproduce the experimentally observed spherical ground states in the 
neutron-deficient Pb isotopes. In RMF, these isotopes are found to be deformed using standard forces with constant 
pairing gaps [10,12]. It should be emphasized that the same RMF-I-BCS calculations that lead to good agreement with 
experimental data in systematic studies of ground state properties carried out over 1300 even-even isotopes [12], fail 
to account for the spherical ground states in this neutron-deficient Pb region. It has been also shown [11] that using 
constant strengths for the pairing interaction, which makes the gap parameters dependent on deformation, produces 
spherical ground states in the Pb isotopes. Yet, more recently [5], an improved treatment of pairing was used by 
means of a relativistic Hartree-Bogoliubov calculation with the NL3 effective interaction and a finite range Gogny force 
to describe the pairing properties. However, it produced again oblate ground states in i88-i94p|^^ contradicting the 
experimental data. In order to recover the spherical ground states, a new parametrization of the effective interaction 
was proposed [5]. 

In a recent work [13], we have calculated the /3-decay properties of the neutron-deficient Pb isotopes using a 
deformed Skyrme HF-I-BCS-I-QRPA approach. This approach was used extensively [14] to demonstrate that the /3- 



decay properties of unstable nuclei may depend on the nuclear shape of the decaying parent nucleus and to predict 
to what extent the GT strengths may be used as fingerprints of the nuclear shapes. Accurate measurements of the 
GT strength distributions in Kr and Sr isotopes [15] have supported the usefulness of these studies. By analyzing 
the energy deformation curves corresponding to those Pb isotopes we found that, although the relative energies of 
the various minima are dependent on both the Skyrmc and pairing forces, the existence of shape isomers as well as 
the location of their equilibrium deformations are rather stable. We also found that the Gamow- Teller (GT) strength 
distributions calculated at the various equilibrium deformations exhibit specific features that can be used as signatures 
of the shape isomers, and what is important, these features basically remain unaltered against changes in the Skyrme 
and pairing forces. Therefore, although an accurate calculation of the excitation energies of the 0"*" states is beyond 
the scope of the theoretical framework used here, our estimate of the /S^-decay patterns is reliable and has been 
proved to be useful [14,15]. 

Motivated by the interest in this type of calculations, not only in Pb isotopes but also in the neighbor regions where 
shape coexistence is found as well, in this paper we extend the calculations of the GT strength distributions to the 
neutron-deficient I80,i82,i84,i86,i88,i90,i92jjg^ I84,i86,i88,i90,i92,i94pb and i98,200,202,204,206po ^^^^ isotopes, using the 

Skyrme force SLy4 [16], which is one of the most recent and successful forces. The reason for this choice of isotopes is 
that they are the best candidates for /5-decay experimental studies. They are /3^-unstable isotopes, yet not dominated 
by other decay modes such as a-decay. In addition they have relative large half-lives and large enough Qec energies 
to make experiments feasible. Nevertheless, we present the results obtained not only within the Qec window relevant 
for the P^ / EC-decay, which is being considered as a real possibility to be measured at ISOLDE/CERN [17], but 
also the GT strength distributions for GT~^ and GT^ in the whole range of excitation energies up to 25 MeV that 
could be extracted from charge exchange reactions in the near future at the new experimental facilities involving 
radioactive isotope beams, such as FAIR. Charge exchange reactions at small momentum transfer are a powerful tool 
to study GT strength distributions when /3-decay is not energetically possible. For incident energies above 100 MeV 
the isovector spin-flip component of the complex effective interaction is dominant and the cross section at forward 
angles is proportional to the GT strength [18]. 

The paper is organized as follows. In Sec. II we briefly present the main features of our theoretical framework. Sec. 
Ill contains our results on the energy deformation curves and Gamow- Teller strength distributions in the neutron- 
deficient Hg, Pb, and Po isotopes relevant for /3+-decay. We also discuss the half-lives, the summed GT strengths, 
and the GT strength distributions in the whole range of excitation energies. Sec. IV contains the main conclusions. 

II. BRIEF DESCRIPTION OF THE FORMALISM 

The theoretical formalism used to describe the GT strength distributions has been already shown elsewhere [14]. 
Here we only summarize the basic ingredients of the method and discuss our choice of interactions and parameters. 
The method starts from a self-consistent deformed Hartrce-Fock mean field calculation with effective two-body density- 
dependent Skyrme interactions including pairing correlations in BCS approximation. In this paper we consider the 
force SLy4 [16], which has been specially designed to improve the isospin properties of nuclei away from beta stability 
and therefore, it may be more adequate to describe the neutron-deficient Hg, Pb, and Po under study in this work. 
Comparison with other widely used Skyrme forces like Sk3 [19] and SG2 [20] was done in the past [13,14] concluding 
that the main characteristics of the /3-decay patterns are not significantly changed when using different forces. 

In this work we have restricted our study to axially deformed nuclear shapes. The single-particle wave functions 
arc then expanded in terms of the eigenstates of an axially symmetric harmonic oscillator in cylindrical coordinates, 
using eleven major shells. This simple picture is validated from triaxial calculations [3,6,8], which show at most slight 
triaxiality (7 ^ 10°) only in very few cases. However, the triaxial barriers are not in general sufficiently high to fully 
support the axial symmetry assumption and it would be interesting to consider these effects in future works. 

The method also includes pairing between like nucleons in the BCS approximation with fixed gap parameters for 
protons A^, and neutrons A^, which are determined phcnomcnologically from the odd-even mass differences through a 
symmetric five term formula involving the experimental binding energies [21]. When the experimental masses are not 
known we use the expression A = 12A^^/^ MeV for protons and neutrons. Occupation probabilities are determined at 
each Hartree-Fock iteration by solving the corresponding gap and number equations. The energy curves as a function 
of deformation are calculated performing constrained HE calculations with quadrupole constraints [22], minimizing 
the HE energy for each nuclear deformation. 

We add to the mean field a separable spin- isospin residual interaction to describe GT transitions. The advantage 
of using separable forces is that the QRPA energy eigenvalue problem is reduced to find the roots of an algebraic 
equation. The residual interaction contains two parts, particle-hole (ph) and particle-particle (pp). The ph part 



y^T = '^xtT E {-^tP\f^--K. ^\^Y.^v\oK\^)ata^, (2.1) 

K=0,±1 -Kv 

is responsible for the position and the general structure of the GT resonance [14,23,24]. The usual procedure to 
determine the coupling strength x^j, is to reproduce the energy of the resonance [24] . In this work we use the same 
value of ')<^Qrj^ for all the isotopes considered in this mass region. The value is fixed by comparison of the calculations in 
^"^Pb to the experimental GT strength resonance observed at an energy of 19.2 MeV, relative to the parent nucleus, 
from (p, n) charge exchange reactions [25]. In this way, we reproduce the energy of the resonance with Xgt ~ ^.08 
MeV. 

The particle-particle residual interaction is a neutron-proton pairing force in the J" = 1+ coupling channel. We 
introduce this interaction as a separable force in the usual way [24,26] 

Vir^ = ~2k^S^ Y,{-1)^P+P^k. Pi^Y. {^\i''K)^\ ^) «>* • (2-2) 

Once the value of the pft, spin-isospin residual force has been established, the coupling strength k^qj- is usually adjusted 
to reproduce the measured half-lives [24,26] because the peak of the GT resonance is almost insensitive to the pp 
force. In our case, as explained in Ref. [13] we have found that the calculated half-lives as a function of kFqj- are very 
little dependent on this parameter and we use k^^, = 0.02 MeV, which is compatible with the parametrization of Ref. 
[24], obtained from a systematic fitting procedure over the nuclear chart. 

For even-even nuclei, the Gamow- Teller strength B{GT^) in the laboratory frame for a transition I^^Ki{0^0) -^ 
1/ Kf{l~^K) can be obtained as 

B{GT^) = Yl K^f^^f |/?^l ^^^^^)r = {^K,fi {(f>K, |/3o^| 0o)' + ^Sk,s {(f>K, |/3i^| 0o)'} , (2.3) 

Mi,Mf,fi 

in units oi g^/Air and in terms of the intrinsic amplitudes calculated in QRPA [14]. To obtain this expression we have 
used the initial and final states in the laboratory frame expressed in terms of the intrinsic states {(pK >, using the 
Bohr-Mottelson factorization [27]. 

The /3-decay half-life is obtained by summing up all the allowed transition probabilities, weighed with phase space 
factors, up to states in the daughter nucleus with excitation energies lying within the Q^c-window, 

^i72 = ^E/(^'^)^(G^)' (2-4) 

where D — 6200 s. We include the standard effective factor 

A^ = [{gA/gvUf = [0-77 {gA/gv),,J ■ (2.5) 

In this work we use experimental Qec values because as it was discussed in Ref. [13] the Qsc-values calculated 
from the binding energies of parent and daughter nuclei are quite similar for the various shapes and are close to the 
corresponding experimental values. In /3+/SC decay, / {Z,lij) consists of two parts, positron emission and electron 
capture. We have computed them numerically for each value of the energy. The Fermi integrals f^ (Z, w) are given 

by 

f^ {Z,Wo)^ pW{WQ-WfX^{Z,W)AW, (2.6) 



with 



A±(Z, W) ^ 2(1 + ^)(2pi?)-2(i-7)eT-v£|l±M! , (2.7) 



where 7 — -^/l — (aZ)^ ; y ~ aZW/p ; a is the fine structure constant and R the nuclear radius. W is the total 
energy of the /3 particle, Wq is the total energy available in TTigC^ units, and p = \/W^ — 1 is the momentum in rrieC 
units. In the numerical calculation, we have included screening and finite size effects as explained in Ref. [28]. The 
electron capture rates f^"-^ have also been included following Ref. [28]: 

X 

where x denotes the atomic subshell from which the electron is captured, q is the neutrino energy, g is the radial 
component of the bound state electron wave function at the nucleus, and B stands for other exchange and overlap 
corrections [28]. 



III. RESULTS AND DISCUSSION 

In this section we present the results obtained for the /3-decay patterns of the neutron-deficient i80-i92jjg^ i84-i94pj-| 
and i98-206pQ gven isotopes. First, we discuss the energy surfaces and shape coexistence expected in these isotopes. 
Then, we present the results obtained for the Gamow- Teller strength distributions with special attention to their 
dependence on the nuclear shape and discuss their relevance as signatures of deformation to be explored experimentally. 
Finally, we discuss the half-lives and the summed GT strengths both within the Qec and in the whole energy range. 

A. Equilibrium deformations 

We show in Fig. 1 the HF energies calculated with the Skyrme force SLy4 and pairing correlations treated in the 
fixed gap approach. We show the results for Hg (left), Pb (middle), and Po (right) isotopes as a function of the 
quadrupole deformation parameter calculated as /3 = ^Jn /bQp/{Z < r^ >), which is defined in terms of the proton 
quadrupole moment Qp and charge r.m.s. radius < r^ >. The curves have been scaled to the energy of their ground 
states and have been shifted by 1 MeV from one isotope to the next one, starting from the lightest one. 

In the case of Hg isotopes we get prolate and oblate minima in all the isotopes from A — 180 up to yl = 192. 
The ground state is prolate for i*°Hg and ^^^Hg and oblate for i84-i92jjg isotopes. Comparable results were also 
obtained in RMF calculations [12], where oblate ground states where obtained at /3 ~ —0.2 with prolate minima 
close in energy at /3 ~ 0.13. This is in qualitative agreement with experiment, where it is observed that a weakly 
oblate ground-state band is crossed in the vicinity oi A = 188 by a deformed band associated with a prolate energy 
minimum. The prolate states minimize their energies for A = 182, although they still lie above the ground state. 
We get qualitative agreement although the details of the energy minima depend much on the effective Skyrme force 
and pairing treatment used in the calculations. In the case of Pb isotopes we obtain a prolate ground state for ^^''Pb 
and oblate ground states for i86-i94p|-, jgotopes. As we can see in the figure, prolate and oblate minima are always 
present from A = 184 to ^ = 190. Spherical local minima can be clearly seen also in A = 184 and shallow profiles 
are found in other cases. For Po isotopes we find oblate ground states in all the isotopes considered i98-206pQ^ -^g 
also find prolate solutions, which are more prominent in A = 200, 202, 204 isotopes. 

In general, we observe that both oblate and prolate deformed minima occur at smaller deformation as the number 
of neutrons increases, approaching the N = 128 closed shell. This is true in Hg, Pb, and Po isotopes and the effect is 
more pronounced in the prolate shapes. The oblate minima are rather stable at /3 = —0.2 in Hg and Pb isotopes and 
at about /? == — 0.1 in Po isotopes. 

We get similar qualitative results for the Hg, Pb, Po isotopes considered when other Skyrme forces are used. More 
specifically, we obtain the same patterns of coexistence with minima at about the same deformations although the 
relative energies may change from one force to another. In the case of Hg and Pb isotopes the correspondence of the 
minima obtained with different forces is almost perfect. In the case of Po isotopes we find the deformation of the 
minima with SG2 and SLy4 almost coincident. Sk3 shows only one oblate minimum that corresponds quite well with 
the oblate minimum obtained with the other forces, but in the prolate region Sk3 shows only a small distortion of the 
curve without developing a minimum. In the case of ^"^Po, Sk3 and SG2 exhibit only a flat minimum centered at the 
spherical shape. 

In general, our results agree with experiment in the sense that we find different coexisting shapes in this mass 
region, but in some cases they are at variance with experiment in the predicted shape of the ground state. This is 
especially true with respect to Pb isotopes. This discrepancy, also found in different theoretical frameworks [5,9], was 
discussed in [13], where we demonstrated on some examples the sensitivity of the energy curves to fine details of the 
two-body interactions, showing that the deformation minima remain at about the same deformation values but their 
relative energies may change considerably. 

B. Gamow- Teller strength distributions 

In this subsection we study the energy distribution of the Gamow- Teller strengths calculated at the equilibrium 
shapes that minimize the energy of the nucleus. 

In Fig. 2 we show the GT strength distributions in (g^/An) units as a function of the excitation energy in the 
daughter nucleus for the Hg isotopes under study. The results correspond to QRPA with the force SLy4, fixed pairing 
gaps obtained from experimental masses and residual separable interactions ph and pp with coupling strengths given 
by Xgt ~ ^-^^ MeV and n^j, = 0.02 MeV, respectively. The left panels contain the discrete strengths for both 
prolate (upward) and oblate (downward) equilibrium deformations. The dashed arrows indicate the experimental 



Qec values of the decay. We have used the same energy scale to appreciate better the differences as we increase the 
number of neutrons. The scale of energies is 6 MeV, enough to include all the Qec energies, whose maximum value is 
Qec = 5.352 MeV in ^®°Hg. The right panels contain the accumulated GT strength for prolate (thick lines) and oblate 
(thin lines) deformations. In this case the energy axis is extended up to the Qec value in order to observe directly 
the allowed energy window and the total GT strength expected in the decay. Although such kind of accumulated 
plots contain in principle the same information as the pure spectrum, it offers a more global and appealing view, 
which is sometimes preferable. One should also notice that in figures 2-7 the quenching factor (see Eq.(2.5)) is not 
included and therefore, one should consider a reduction of this strength by a factor of about two when comparing to 
experiment. 

From Fig. 2 one can see that prolate and oblate nuclear shapes give rise to GT strength distributions whose 
patterns can be clearly distinguished from each other. The oblate strength appears in all cases very fragmented over 
the whole Q-window, which induces a straight increasing of the strength in the accumulated plot. On the other hand, 
the distributions in the prolate cases show a concentration of the GT strength in a single strong peak located at very 
low excitation energy (below 0.7 MeV in all the isotopes). This single excitation carries practically the whole strength 
contained in the Q-window in the A = 190, 192 isotopes and about half of it in ^ = 180 — 188. The accumulated 
plots are then characterized by the successive appearance of large vertical steps followed by large flat regions, in 
contrast to the oblate case. With the help of the accumulated plots we can also observe immediately that although 
the total strength contained within the Qec energy is quite similar for oblate and prolate shapes, the distribution of 
this strength along the energy axis is very different. 

Fig. 3 contains the same information as Fig. 2 but for Pb isotopes. In this case we also include results for 
spherical shapes in such a way that the left panels contain the results from prolate shapes (upper plots), oblate shapes 
(middle plots) and spherical shapes (lower plots), except in the A ~ 194 isotope where we only have two shapes. 
Similarly, the right panels contain three lines, thick line for prolate, thin for oblate, and dashed for spherical shapes. 
As in the case of Hg isotopes, the signature of the prolate shapes is a strong isolated peak at low excitation energies 
(there is a second strong peak in A = 188). This makes again the accumulated plots to appear as big steps. The 
oblate distributions show a very fragmented structure, starting at higher energies than the strong prolate peak. This 
produces an accumulated pattern which increases in a steady way, starting at higher energy than the prolate pattern. 
In the spherical case, the strength is practically collected in a single peak, which corresponds to the 7r/iii/2 -^ ^^^^9/2 
transition. This excitation carries less strength than the prolate case and occurs in all the isotopes except in A — 184 
at an energy higher than the prolate peaks. It is also worth noticing that oblate and prolate shapes generate practically 
the same total GT strength within the Qec window although the internal distribution of this strength is in general 
very different. The total strength in the spherical case is always clearly below the prolate or oblate strength, except 
in the isotope A = 194, where the small strength and Qec rnakes it not very relevant. 

Fig. 4 shows the same results as in Fig. 2 for Po isotopes. In the cases oi A — 198, 200 we find in the prolate 
case below Qec a strong single peak carrying most of the strength, while again the oblate shapes produce a very 
fragmented distribution. In the other cases, A — 202,206, the strength below Qec is insignificant and could not tell 
much to distinguish the shapes. 

In general, we observe that the GT strength distributions corresponding to prolate shapes in both Hg and Pb 
isotopes are characterized by the appearance of a strong excitation below 1 MeV that collects a large amount of the 
total strength contained in the whole region below the Qec energy. On the other hand, the GT strength distributions 
corresponding to oblate shapes are always very fragmented and are extended over the whole Qec window. In the case 
of the Pb isotopes this strength is located at excitation energies higher than the prolate peak. For Pb isotopes the GT 
strength distributions corresponding to spherical shapes appear as single peaks carrying all the strength below Qec 
and also at higher energies than the prolate peak. These features are in agreement with the conclusions in Ref. [13] 
using the force Sk3. The most noticeable difference with respect to the results obtained from Sk3 force [13] appears 
in the spherical case, where the GT strength obtained with Sk3 is displaced to higher energies as compared to the 
strength obtained with SLy4. The origin of the shift in the excitation energy can be traced back to the predicted 
energies for the ft-g/2 and /111/2 spherical shells. The main contribution to the GT strength at low energies comes 
from the transition connecting the almost fully occupied T^hiii2 shell with the partially unoccupied 1^/19/2 shell, and 
the relative position of the two shells determines the GT excitation energy. The effect of deformation is to allow for 
multiple transitions fragmenting the strength and smoothing the differences caused by the forces. 

In Figs. 5-7 we show the B{GT^) strength distributions in the Hg, Pb and Po isotopes. We plot continuous 
distributions resulting from a folding procedure using F = 1 MeV width gaussians on the discrete spectrum. The 
energy range of excitation energies in the daughter nucleus is extended up to 20 MeV in order to cover all the strength. 
We find the strength beyond 20 MeV to be insignificant, except for Po isotopes where we extend the scale up to 25 
MeV. The interest of these figures, not limited by the Qec energies, is to show the whole strength distributions that 
could be explored from charge exchange reactions such as (n,p) and {p, n) or similar reactions. The B{GT^) strength 
on the right panels is much stronger than the corresponding B(GT^) strength on the left panels, as it should be 



according to the Ikeda sum rule T,^[B{GT~){uj) — B{GT^){uj)] — 3{N — Z), which is fulfilled in all our calculations 
up to a small few percent discrepancy. Ikeda sum rule ranges from 3{N — Z) — 60 in ^*°Hg up to 3(iV — Z) — 114 in 
206pQ^ We can see in the last three columns of Tables 1, 2, and 3 the total sums of the B(GT^), B{GT^), and the 
percentage of the Ikeda sum rule fulfilled in the calculation for Hg, Pb, and Po isotopes, respectively. 

In Fig. 5, the profiles of the B(GT^) strength show a first bump, which for the lightest isotopes considered is a 
double bump, moving from 4 McV in A = 180 up to 7 McV in A = 192. We can also see a big resonance at higher 
energies, that moves from 11 MeV in A — 180 up to 14 MeV in A ~ 192. In any case, the profiles obtained from oblate 
and prolate shapes are quite similar and then, it will be very difficult to use this information to distinguish between 
one shape or another. In the case of the B{GT'^) strength, there is a resonance at low energy below 2 MeV, which is 
rather constant in all the isotopes. As we have seen in more detail in Fig. 2, it is in this region below Qec where the 
shape dependence could be exploited because prolate and oblate shapes produce distinguishable distributions of the 
GT strength. 

Fig. 6 is similar to Fig. 5 for Pb isotopes. Again we see in the profiles of the B{GT^) strength on the right 
panels a small bump between 5 and 6 MeV and a resonance centered at 12-13 MeV. The nuclear shapes can hardly 
be identified looking at the strength distributions. The B{GT'^) strengths show also a resonance at low energies, 
which corresponds to the analysis made in Fig. 3. Fig. 7 for Po isotopes shows once more the small bump in the 
B(GT^) distribution at 6-7 MeV and the big resonance at 12-14 MeV with similar characteristics in all the isotopes 
and shapes. The B{GT'^) strength shows in this case a resonance at high energy around 20 MeV and the low energy 
distributions that were studied in Fig. 4. 

C. Half-lives and summed strengths 

It is also interesting to calculate integral magnitudes from the GT strength distributions that characterize them in 
a single quantity. This is the case of the total half- lives and total GT strength contained in the energy range available 
in the /3+-decay. These magnitudes depend differently on the GT strength distribution, one of them is simply a sum 
and the other one is a weighted sum (see Eq.(2.4)). While the total GT strength involved in the /J^-decays are not 
measured yet, there is experimental information on the half-lives [21] that allows to contrast our calculations and to 
check that there is no substantial disagreement with experiment. We can see the results in Tables 1, 2, and 3 for Hg, 
Pb, and Po isotopes, respectively. The first column in Tables 1-3 contains the experimental Qec energies [21]. Since 
the f3~^ / EG decay mode competes in this mass region with a-decay, we show in the second column of the tables the 
percentage of the total decay assigned to /3+ [21]. The third column contains the total experimental half-lives [21] 
and within brackets the corresponding /J"*" half-lives extracted from the percentage in the second column. The next 
columns contain theoretical results calculated at the equilibrium shapes obtained with the force SLy4. We can see 
first the (3^ /EC half-lives calculated with quenching factors (2.5) and experimental Qec values, as discussed in Ref. 
[13]. Next, we can see the B{GT'^) strength summed up to the Qec energies, the total B{GT^) and B{GT^) sums 
and the percentage of the Ikeda sum rule fulfilled. 

In Table 1 one can see that the half-lives obtained for Hg isotopes agree in general with experiment within a factor 
of 2. Only the most stable isotope A — 1%2 deviates considerably from this agreement, but this case is not very 
relevant because of the small energy window allowed for the decay Qec = 0.765 MeV and because of the small 
strength involved, which makes it meaningless. A factor of 2 is perfectly acceptable in this type of calculations taking 
into account the uncertainties coming from the forces used, the level of approximations, the quenching factors and 
the Qec energies. For instance, reasonable quenching factors could change the half-lives within a 25% effect and 
varying Qec values within 1 MeV range could change the half-lives by a factor of two. Another source of uncertainty 
in the half-lives is induced by the determination of the pairing gaps from the experimental binding energies. Through 
variations of 20% below and above the adopted proton and neutron pairing gap values, we obtain half-lives which in 
general do not differ by more than 20% from the original results. We also find that in the spherical cases the half-lives 
are more sensitive to these uncertainties than in the deformed cases. Nevertheless, as we have discussed, this does not 
rule out the conclusions extracted based on the patterns shown by the GT strength distributions. We also include 
in the table the strength contained below the Qec energy that could be measured in /?-decay experiments [17]. We 
should keep in mind that in contrast to the calculation of the half-lives, we do not include quenching factors in the 
summed strengths and thus, a reduction of about 50% of the values written in the tables are expected when comparing 
to experiment. The sums over the whole energy region of the B(GT^) and B{GT~) strengths show that the total 
B{GT~) is much stronger than B{GT'^) as it should happen in order to fulfill Ikeda sum rule, as it was discussed 
above. In the case of Hg isotopes we can see that the half-lives from prolate deformations are in general larger than 
those from oblate ones with the already mentioned exception of A = 192. The summed strengths are practically the 
same no matter what the deformation is. 



The same general comments can be applied also to Table 2 for Pb isotopes. We get in general half-lives that agree 
with experiment within a factor 2-3 . This is also true when we consider other Skyrme forces. In particular, the 
half-lives obtained with the forces Sk3 and SG2 are very similar to the corresponding SLy4 ones shown in the tables. 
The only disagreement appears in the spherical cases. We find in the spherical case differences up to factors of 4-5 in 
the half-lives obtained from various forces. As we have already discussed above, this is due to the different predicted 
energies of the spherical shell gaps. If we compare the half-lives from different deformations, we can sec that the half- 
lives from the spherical shapes are larger than those corresponding to the deformed shapes, and the prolate shapes 
give always lower half-lives than the oblate ones. In the case of SLy4, the half-lives from the spherical shapes agree 
with experiment at the same qualitative level as the half-lives from deformed shapes. This is at variance with what 
we found in Ref. [13] using Sk3 and SG2 forces, in the sense that there, we only found agreement with experimental 
half-lives with deformed shapes. The strengths summed up to Qec E^re very close to each other for oblate and prolate 
deformations and arc much larger than the corresponding strength for the spherical case. 

In Table 3 we have the results for Po isotopes. In this case the total GT strength involved in the decay is very 
small in all the cases and therefore the half-lives are large and not very meaningful. Only in the two more unstable 
isotopes A = 198, 200 this comparison makes sense and we can see that again the half-lives agree with experiment 
within a factor of two. The oblate half-lives are larger than the corresponding prolate ones. 

Although the deformation does not show up in the total GT strengths, it docs in the internal distribution of the 
GT strength, specially in the B{GT^) at low excitation energy where the /3+-decay is sensitive. The only purpose of 
showing half-lives and summed strengths is to contrast the validity of the calculation against some measured quantities 
in the case of half-lives or to predict the expected total strength coming into play in the other case. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the energy distribution of the Gamow- Teller strength in the neutron-deficient i80-i92jjg^ i84-i94p^^ 
and i98-206pQ isotopes. We use a deformed pnQRPA formalism with spin-isospin ph and pp separable residual 
interactions. The ph coupling strength is fixed to reproduce the energy of the GT resonance in ^"^^Pb. The pp 
coupling strength is taken from existing global parametrizations, although we find our results to be quite insensitive 
to this value. The quasiparticle mean field is generated from a deformed HF approach with the two-body Skyrme 
effective interaction SLy4. It includes pairing correlations in BCS approximation, using fixed gap parameters extracted 
from the experimental masses. The equilibrium deformation is derived sclf-consistently within the HF procedure. 

An analysis of the energy deformation curves shows that oblate and prolate coexisting shapes are expected in the 
neutron-deficient Hg and Po isotopes. In the case of Pb isotopes we also get in some cases a spherical minimum 
giving rise to a triple shape coexistence. We should remark again that the relative energies of the minima in these 
calculations are very sensitive to the Skyrme force and type of pairing used and therefore the predicted shapes for the 
ground states can be altered depending on the choice of these forces. However, for a given isotope, the equilibrium 
deformations at which the minima occur are hardly shifted. 

The GT strength distributions calculated at the equilibrium shapes show that the effect of the deformation is 
much stronger than the effects coming from the Skyrme or pairing forces used. As a general rule, we find that the 
GT strength distributions calculated from prolate shapes are characterized by strong single peaks at low excitation 
energies carrying a sizable amount of the total GT strength involved in the /3^-decay. The small remaining GT 
strength is scattered over the energy interval up to Qec- On the other hand, the signature of oblate shapes is a 
distribution spread all over the available Q- window energy. In the case of Pb isotopes, the GT strength distributions 
calculated from spherical shapes are characterized by single peaks that contain all the GT strength, which are in 
general located at higher energies than the prolate peaks. The detailed analysis performed allows us to identify good 
candidates to explore experimentally the signatures of deformation on the GT strength distributions. We find that 
the most interesting isotopes are i84-i88jjg^ i88-i92p^ ^^^^ i98,200py^ 

We have also investigated the GT strength distributions in the whole range of excitation energies without the 
limitations imposed by the /3"'"-decay. This analysis is worth as a prediction for the cross sections to be expected in 
charge exchange reactions. 

The detailed analysis at low excitation energies in the B(GT^) distribution shows that the /3+-decay process is the 
most sensitive to the nuclear shape. 

The calculated total strength involved in the /3~''-decay shows little dependence on the nuclear shape, except in the 
case of the spherical shapes in Pb isotopes, where this strength is much smaller. The calculated (3^ /EC half-lives 
agree with experiment within a factor of 2 or 3. In general we can conclude that neither the half- lives nor the summed 
strengths are good observables to study deformation effects, but the strength distributions are. 
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TABLE I. Half-lives and B{GT ) strengths in Hg isotopes. The table contains experimental Qec values [MeV], percentage 
of the P'^ /EC involved in the total decay, total experimental half-lives [s] [21] and within brackets, the j3'^ /EC experimental 
half-lives extracted from the percentage. Then we find theoretical results obtained with the force SLy4: half-lives [s], B{CT^) 
strength [pA/(47r)] summed up to Qec energies, total B(CT'^) and B{CT~) strengths contained in the whole energy range 
considered {Eex < 30) MeV, and percentage of the Ikeda sum rule fulfilled in our calculations. 



isotope 



lEC 



%I3^ 



Ti/irLp (/3+) 



t: 



■I3+/EC 



1/2, SLyi 



,B{GT+) T,B{GT+) 



T,B{GT- 



% Ikoda 



i«OHg, obi 
prol 

i^^Hg, obi 
prol 

is-^Hg, obi 
prol 

i^'^Hg, obi 
prol 



prol 

i^°Hg, obi 
prol 



5.352 



4.725 



3.970 



3.176 



^Hg, obi 2.099 



1.511 



192 



Hg, obi 0.765 
prol 



52 



86.2 



98.89 



100 



100 



100 



100 



2.6 (4.9) 


3.0 


3.10 


4.63 


64.20 


99.28 




4.6 


2.73 


4.43 


63.99 


99.27 


10.8 (12.6) 


5.6 


2.49 


4.27 


69.62 


99.02 




7.5 


2.28 


4.18 


69.51 


98.98 


30.6 (30.9) 


17.5 


1.71 


3.06 


74.46 


99.17 




17.0 


1.46 


3.08 


74.50 


99.19 


82.8 (82.8) 


47.2 


1.22 


2.57 


79.85 


99.08 




68.4 


0.89 


2.46 


79.76 


99.10 


195 (195) 


185.1 


0.78 


2.51 


85.64 


98.96 




218.9 


0.74 


2.59 


85.72 


98.96 


1200 (1200) 


506.9 


0.63 


2.10 


90.84 


98.60 




493.8 


0.47 


2.27 


91.35 


98.98 


7460 (17460) 


5200 


0.14 


1.71 


96.55 


98.79 




1970 


0.27 


1.80 


96.94 


99.10 
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TABLE II. Same as in Table 1 for Pb isotopes. 



isotope Qec % W^/EC) 



p/S+Z-BCn 



^l/2,exp \-'-l/2,exp 



Ti 



't3+/EC 



l/2,SLy4 



^qec'B{gt+) t.b{gt+) j:b{gt- 



^''Pb, sph 


5.84 


obi 




prol 




*«Pb, sph 


5.51 


obi 




prol 




**Pb, sph 


4.53 


obi 




prol 




''Opb, sph 


3.92 


obi 




prol 




92pb, sph 


3.32 


obi 




prol 




''■'Pb, sph 


2.62 



20 



90.7 



99.6 



100 



100 



0.49 (2.45) 



4.82 (8.03) 



25.5 (28.1) 



71 (71) 



210 (210) 



720 (720) 



obi 



5.6 


0.85 


1.95 


61.65 


99.50 


5.6 


2.31 


3.61 


63.31 


99.50 


3.3 


2.55 


4.01 


63.70 


99.48 


16.3 


0.54 


1.65 


67.28 


99.44 


10.7 


1.73 


2.91 


68.58 


99.50 


4.8 


1.91 


3.34 


68.94 


99.39 


28.0 


0.76 


2.11 


73.76 


99.51 


20.3 


1.35 


3.05 


74.72 


99.54 


16.1 


1.46 


3.16 


74.78 


99.47 


80.0 


0.54 


1.92 


79.48 


99.44 


41.7 


0.99 


2.67 


80.26 


99.47 


26.4 


0.99 


2.89 


80.32 


99.27 


251.7 


0.41 


1.90 


85.57 


99.61 


93.4 


0.70 


2.44 


86.12 


99.62 


45.6 


0.79 


2.70 


86.26 


99.48 


655.8 


0.62 


2.02 


90.70 


98.53 


482.5 


0.56 


2.32 


90.99 


98.52 



TABLE III. Same as in Table 1 for Po isotopes. 



isotope Qec % (P^/EC) 



rptotal frjil3+/EC\ 

-'-l/2,exp \-'-l/2,exp ) 



T, 



11+ /EC 
1/2, SLyA 



-•Qec 



B(GT^ 



T,B{GT^ 



T,B{GT- 



^^**Po, obi 3.900 
prol 

2°°Po, obi 3.416 
prol 



202 



Po, obl 2.809 
prol 



^o-'Po, obl 2.334 
prol 



206 



Po, obl 1.846 
prol 



43 



99.34 



94.55 



106 (247) 


342 


0.387 


2.23 


91.13 


98.78 




160 


0.333 


2.24 


91.25 


98.90 


690 (776) 


1390 


0.085 


1.82 


96.64 


98.77 




1066 


0.189 


1.82 


96.76 


98.90 


2682 (2735) 


8799 


0.029 


1.42 


102.9 


99.46 




5341 


0.015 


1.41 


102.88 


99.48 


12708 (12729) 


29788 


0.013 


1.32 


108.88 


99.59 




12167 


0.018 


1.31 


108.90 


99.62 


760320 (804145) 


256650 


0.0037 


1.25 


114.78 


99.59 




55352 


0.0046 


1.24 


114.72 


99.54 
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FIG. 1. HF energy scaled to the ground state energy (see text) obtained from constrained HF+BCS calculations with the 
Skyrme force SLy4 and fixed pairing gap parameters as a function of the quadrupole deformation (3. Left: i80-i92-pjg jgQ^Qpgg^ 



middle: i«4-i94p,^ isotopes, and right: i98-206po isotopes. 
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FIG. 2. Left: Gamow- Teller strength distributions [5^/(4"")] in Hg isotopes for prolate (upward) and oblate (downward) 
shapes. The experimental Qec energies are shown by dashed arrows. Right; Accumulated Gamow- Teller strength for prolate 
(thick lines) and oblate (thin lines) shapes plotted up to the Qec energies. Results are obtained from SLy4 force with fixed 
gap parameters. 
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FIG. 3. Same as in Fig. 2 for Pb isotopes. In this case we also include results from spherical shapes. 
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FIG. 4. Same as in Fig. 2 for Po isotopes. 
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FIG. 5. Folded Gamow- Teller strength distributions B{GT ) [pA/(47r)] for the various shapes of the Hg isotopes in the range 
of excitation energies up to 20 MeV. 
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FIG. 6. Same as in Fig. 5 for Pb isotopes. 
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FIG. 7. Same as in Fig. 5 for Po isotopes. 
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